PENALTY APPROXIMATION FOR NON SMOOTH 
CONSTRAINTS IN VIBROIMPACT 

LjETITIA PAOLI and MICHELLE SCHATZMAN 

Abstract. We examine the penalty approximation of the free motion of a 
material point in an angular domain; we choose an over-damped penalty ap- 
proximation, and we prove that if the first impact point is not at the vertex, 
then, the limit of the approximation exists and is described by Moreau's rule 
for anelastic impacts. 

The proofs rely on validated asymptotics and use some classical tools of 
the theory of dynamical systems. 



1. Introduction 

Mathematical results relative to the convergence of a penalty approximation 
of impact problems have been obtained by several authors when the energy is 
conserved; see for instance ||, and also i,§i,[|l,§,i,0,i, [| 
and 0. 

When energy may be lost at impact, the convergence of the penalty approx- 
imation has been treated in in the case of a convex set of constraints with 
smooth boundary. In this article, we defined a penalty approximation for which 
the limit solution satisfies a Newton condition at impact: the normal component of 
the velocity is reversed and multiplied by a restitution coefficient e g]0, 1] and the 
tangential component is transmitted. 

So far, we are not aware of any mathematical results on the convergence of the 
penalty approximation when the boundary is not smooth and energy can be lost 
at impact. 

Here, we study the penalty approximation of the motion of a free particle con- 
strained to stay inside an angular domain of R^: we choose a class of penalty 
approximations for which the restitution coefficient vanishes in the limiting prob- 
lem and we characterize precisely the limit of the sequence of solutions of the 
approximated problem when the first impact does not take place at the corner. 

We compare our results to the ones given by the selection rule of Moreau 
and we find complete agreement. 

2. The first part of the motion and the mathematical strategy 

Let us describe more precisely the problem and the method of solution. 
Given 6 G (0, tt), we let K be the set 

(2.1) K = {{xi,X2) e : xi < and xi cos^-h X2 sin^ < O}. 
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Figure 1. The regions and the coordinates in the plane. 



The closure of the complement of K is partitioned into three regions: 

TZi — {{xi,X2) ■ xi >0 and X2 < 0}, 
T^2 — {{xi,X2) ■ X2 >0 and — xi sin^ + X2 cos9 < 0} 
7^3 — {{xi,X2) : xi cos9 + X2 sin^ > 
and — xi sin 6 + X2 cos ^ > 0}. 

In each of these regions, the projection onto K , which is known to be a contrac- 
tion, takes different forms: 

'(0,2:2)'^, ifa;e7^l, 
Pkx^{Q, ifa;G7^2, 
{—xi sin 9 + X2 cos 9){— sin 9, cos 6)^ , if x G TZ^. 



The penalty approximation used in 1 11 is defined as follows: we define a function 
G of two arguments it e IR-^ and f e M~ by 

G[u,v) = i \u~ Pku\ 

[0, otherwise. 

Then, the penalized approximation to the impact problem in K , in the absence 
of exterior forces is given by 

(2.2) Uk + 2aVkG{uk,Uk) + k{uk - Pkuu) = 0. 

In this formulation, the number k is the stiffness of the spring which describes 
the reaction of the wall, and the choice of the scale ^/k is the exact choice which 
ensures convergence in the smooth k tends to infinity. Here, we choose 

a > 1: it is the over-damped choice and it will lead to a vanishing restitution 
coefficient as we shall see. 
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Take initial conditions given by 
^2 3^ Uk{0) = iO,X2{0)f, X2(0)<0, 

Ufc(O) = (ii(0),i2(0))^, ii(0) >0,i2(0) >0. 

These initial conditions mean that at time t — 0, the representative point of the 
system is on the boundary of K, in region TZi, and outgoing as well as taking the 
direction of the corner. In particular, this choice of initial conditions means that 
the first impact time is t = 0. 

The roots of the characteristic equation 

+ + 1 = 

of the over-damped equation 

(2.4) y + 2ay + y = 
are given by the formulas 

(2.5) A^a^~l, ^i=-a + VA, ^2 = -a - VA. 

Both ^1 and ^2 are negative. 

As long as the representative point of the system lies in TZi, we perform the 
change of variables r{t) — xi{t) > and s{t) = X2{t) < 0. In these new coordinates, 
( p^ ) becomes the decoupled system 

r + 2aVkr + kr — 0, 
s = 0. 



Its solution is given explicitly by 

(2.6) r{t) - fe«^* 

2VAfc V 

(2.7) s{t) ^ s{0) + ts{Q) 



2VAfc 



For all positive i, r{t) given by ( p.q ) remains strictly positive; s reaches the value 
at the time 

(2.8) to = -mim- 

Therefore, at the boundary between regions TZ\ and 1^2 we have 
s(to-0) = 0, s(to-O) = s(0), 

(2.9) 2VAfcV^ r 
f{to - 0) = ^ (ae«^*«^ - Cae^^*"^) . 

In order to study the motion in region 'R,2^ we use polar coordinates i.e. Uk = re^^ 
and we define scaled functions and variables R, O and t by 

r(t) =77i?(T)/%/fc, T={t-to)Vk, ,y = e^i'o^/^, Q{T)^e{t). 

We have represented in Fig. |^ the numerically computed trajectories (dotted or 
dashed lines) and the vector field of the ordinary differential equation for R and R. 
In the new variables, the system under consideration becomes 

(2.10) R ^ „^ ' +2aR + R = 0. 

R 




with 



(2.11) 



e = 



i?2 



and the detailed derivation of these equations is performed in subsection 3.1. In 
equations (2.1C) and (2.11), e — o(l) and is a fixed number depending only on the 
initial conditions and a. The representation given in Fig.0 will help us to explain 



how the solution of (2.1C) behaves, with appropriate consequences on the angle Q. 

In region Ai, R decreases somewhat and then increases, R increases from a size 
equivalent to Crj to a size equivalent to C /rj in that same region. The dominant 
terms in equation ( 2.10| ) are R and E{1 — e)'^/R^; therefore, we are led to the 
problem 

E 



(2.12) 



Ri 



Rl 



0, i?i(0) = i?(0), i?i(0) = i?(0). 



We study the solution Ri of ( 2.12| ) in subsection |3.2| , as well as the evolution of the 
function Oi satisfying 



this can be done explicitly, thanks to the simple structure of ( ^.12 ). In subsection 
|3.3| , we study the kernel of the linearized (2.12) at as a preparation for the 
validation of this first asymptotic, a task which is completed in 3.4 on the interval 
[0,Ti], where ti is equal to ry'''^, with 71 belonging to (1,2). We conclude section 
H by Proposition 3^ which shows that R is equivalent to Ri over [0,ri] and R is 
equivalent to Ri over [77'^, ri]. The proof is basically a consequence of the fixed 
point theorem with a number of technical estimates. 
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In Section ^, assuming 6 < 7r/2, we are able to exploit the above equivalents 



and to prove that Q, solution of ( |2.1lD , crosses through 9 at some time f < 0[Tf). 
Moreover, our estimates enable us to describe the limit Uqo of Uk as k tends to 
infinity. Let Hi be the orthogonal projection on {xi — 0}, and let 112 be the 
orthogonal projection on {xi cos9 + X2 sui9 = 0}; then 



u(0) + tni?i(0) ifO<<<io, 

(i - io)n2niM(io) if<o<t. 



If 6* > 7r/2, the representative point of the system enters region A2 of Fig. |2[ We 
have to produce an asymptotic for the solution of ( 2.10| ); in this region, it is the 



linear part of this ordinary differential equation which is dominant; more precisely, 
let i?2 be the solution of 

i?2 + 2ai?2 + R2 = 0, 

with i?2 and R2 respectively coinciding with R and R at time ri = rj'^^ , where, now 
the interval of 71 is reduced to (1,4/3). 

The validation of this ansatz is another consequence of the fixed point theorem 
for strict contraction, together with a number of technical estimates. 

Finally, we use classical methods for dynamical systems and prove that the 
representative point of the system tends to {Rc,0) as time tends to infinity: Rc 
is a number which depends only on the initial conditions, a and e. We combine 
the use of a Lyapunov functional and some elementary properties of the system 
to conclude that R remains bounded from above and away from for all time 
after leaving A2. Observe that the Lyapunov functional gives scant information in 
regions A2 and Ai: there it takes values of order l/rj. 

With some technicalities in the case 9 = 7r/2, it is possible to conclude that Q{t) 
crosses 9 at some time f and to obtain precise equivalents for R, R and O at time 



After this time, the representative point of the system (2.2) enters region TZ^, 



and we conclude by Theorem [7^ that the limit Mqo of Uk is given by 



u(0) + <niu(0) ifO<t<io, 
if to < t- 



Moreau's rule is described as follows: at impact, the outgoing velocity is pro- 
jected onto the tangent cone to the convex of constraints, and the motion proceeds 
with this new velocity. Thus, it can be seen that the over-damped penalty approxi- 
mation agrees completely with Moreau's rule if the first impact does not take place 
at the corner, or very close to it, i.e. at a distance 0(l/\/fc) from it. 

We conjecture that the behavior described here still holds if there is a right hand 
side, and the convex is replaced by a set with convex corners, and smooth and not 
necessarily convex curves between corners. We also conjecture that the behavior 
of the limit of the over-damped penalized solution is the same in higher spatial 
dimension. 



3. Equations of the motion around the corner: the earliest 

asymptotic 

3.1. Derivation of the scaled equation in Tl2- After time to, we arrive into re- 
gion TZ2, in which it is convenient to identify and C and to use polar coordinates. 
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ie 



(3.1) Uk = re 

By continuity, the limits of r(t) from the right and from the left as t tends to to are 
identical; therefore: 

(3.2) ^(to + o) = ^^(e«^*"^-e«^*°^). 
At t = to, 

(3.3) e{h+Q) = Q. 
We differentiate once ( ^.l[ ) with respect to time: 



hence 



ilk = re^ + idre 

(3.4) f(to + 0) = ^ (eie«i*°^ _ ^^e^2toVk\ 

(3.5) e{to + o) 



2%/A 

s{to - 0) _ s{0)2VAk 



r{to - 0) ^(-Q) |^e«i*o^ - e^2*oVI^ 
Let us derive the differential equations satisfied by r and 6. In region Tl2, the 



definition of the projection Pk implies that (2.2) can be written as 

(3.6) r9 + 2re = 0, 

(3.7) f - r9^ + 2ay/kr + kr ^ 0. 

The motion has central acceleration, therefore it has a first integral: there exists a 
constant T such that 

(3.8) (r2^)(t) = (r2^)(to + 0) = r. 



and, according to (pA) and (3.5), the value of F is given by 



We substitute 6 = T/r into (3.7), and we find the equation in 



(3.10) r J + 2aVkr + kr = 0. 

Let us find now appropriate changes of scale which will help us understand the 
behavior of r while Uk remains in region TZ2 ■ An important scale is defined by the 
number 

(3.11) ,y = e«i*«^/2^ 

which is very small since is strictly negative. We perform the following changes 
of variables: 

(3.12) T = {t-to)Vk, r{t) =rjR{T)/Vk. 
In these new variables, we have 

f{t) = r]R{T), r{t) = ■qVkRir), 
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SO that equation (3.10) becomes, after a division by rj^/k 
(3.13) 

Let us define 



fcr^ 

R — + 2aR + i? = 0. 



(3.14) 
(3.15) 



E = 



f(0)'s(0)^ 
4A ' 

M2-il)to^/k 



With these notations, we find that F, given by (3^), is equal to 

^ r;2(l-e)v^ 



(3.16) 



Vk 



so that ( 3.13 ) can be rewritten as 
(3.17) 



R ^ ' + 2aR + R = 0. 

R'^ 



The initial values for ( p. 17 ) are given by 

77f(0)(l-e) 



(3.18) 
and 
(3.19) 



i?(0) 



2^/A ' 

77r(0)a(l-£6/a) 



2VA 



(3.20) 



3.2. Ansatz for the motion in Tl2- We shall u se no w an ansatz, namely, we 
state that the essential part in the left hand side of ( 3.17 ) is R — E/ R^ . This comes 
from the fact that at time to, E/R^ is very large with respect to R and R as can 
be checked from (3.18) and (3.19). Therefore, we first solve explicitly the equation 

E 

r! 

If we multiply ( 3.20|) by i?i and integrate, we find that 
(3.21) ^ -^'-^ ^ 

Let us denote by 
(3.22) 



^1-153=0' ^i(0)-i?(0), i?i(0) = i?(0). 



the value which appears on the right hand side of ( ^3.21 ) . Thanks to the change of 
variable p = Rf, equation ( 3.21 ) becomes: 

Wp 



(3.23) 



2y/Wp - E 



±W. 



At the initial time, p{0) is strictly negative, so that in (3.23) we choose the minus 
sign on the right hand side, we integrate until a time tq such that p vanishes, and 
we find 



(3.24) 



E + {Wr + R{0)RiO))^ 

P{T) = , < T < To. 
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The value of tq is given by 
(3.25) 



To = -R{0)R{0)/W. 



After Tq, we choose the plus sign in (3.23), and we find that 
(3.26) p{t)W - E ^ (t - To)'^W'^ . 



Substituting the value of tq into (3.26), we find that the general expression of the 
solution of (3.23) is given by 



(3.27) 



p{r) = § + W{r-Tor. 



The angle 9 is defined by (^); we let 

(3.28) e{T)^9{t). 

We are only interested for the present moment in the principal part of 0; it is a 
function Oi which satisfies the ordinary differential equation: 



(3.29) 

with the initial condition 
(3.30) 



ei(r) = ^, 

p 



ei(o) = 0. 



We substitute the value of p given by (3.27) into (3.29) and we find that 



(3.31) 



ei(r) 



E + {Wt + R{0)R{0)y M-r? 

which we integrate immediately into 

W^T + -R(0)i?(0) 7?(0)i?(0) 
ei(r) arctan _ ^ ^ ' - arctan - ^ ' ^ ' 



E 



E 



A more convenient way to write Oi is the following 
(3.32) 



9.(r) ^ arctan arctan^-" 



E 



E 



Thanks to ( |3.18| ) and (|t|), we can see that 

(3.33) W^s{Qfri-'^. 
But ( p33| ) and ( |3^ ) imply that 

(3.34) To = 0(7,4). 

This shows that we shall have to consider different cases: 9 < tt/2 and 6 > tt/2. 
If ^ < vr/2, we may suspect that Uk will exit region TZ2 at time approximately 
to + [-s/E tan (Pj / (W Vk) ; while if ^ > n/2, it is obvious that the ansatz is not 
sufficient: what will happen is that E/R^ is no more large with respect to R and 

it. 
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3.3. Kernel for the linearized equation. In order to go further, we have to 
vahdate asymptotics: consider therefore the hnear differential equation obtained 
from the linearization of (3.2C): 

(3.35) ^+^ = 0. 



2,Ez 



We can obtain two linearly independent solutio ns by the following argument: if we 
differentiate the ordinary differential equation ( 3.20 ) with respect to time, wc find 
that zi = 7?i is a solution of ( 3.35| ); we seek another solution of ( |3.3E| ) under the 
form 



(3.36) 



Z2 = ZiT. 



indeed, formula (3.27) gives the following form of Ri 



E 



(3.37) 

We differentiate this relation with respect to r, and wc find that 
(3.38) i?i(r)=^("-"°) 
The equation satisfied by T is 



R^ 



2ziT + ziT = 0. 

We multiply this equation by zi, we integrate, and we find that, up to an irrelevant 
multiplicative constant, T satisfies the ordinary differential equation: 

1 

This can be integrated exactly and we obtain 

E T 
V3(T-ro) ^ W' 



T ■ 



T 



According to definition (3.36), we have to multiply the above expression by Ri, for 
which we take expression ( 3.381 ). We obtain eventually 

E r(T-To) 
W^Ri Ri ■ 

The Wronskian of zi and Z2 is readily calculated and is equal to 



^12:2 



Z2Z1 = zlt = 1; 



therefore zi and Z2 are independent. From here, we seek a kernel K{t^ a) for t > a 
which satisfies the following conditions: 

' d^K{T, a) 3EK{t, a) 



9t2 Rf{T) 

K(a,a) =0, 
dK{T,a) I 



= 0, 



dr 



1, 



under the form 



K{T,a) ^ ai{a)SiiT) + a2{<j)S2{T). 
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Thanks to ( |3.3| ) and the definition of zi, we can see now that 

K{t, cr) = -Z2((7)zi(t) + Zi(ct)z2(t), 

which can be rewritten as 



(3.39) 



K{t, a) 



E 
W 



Ri{T)Ri{a) 

We extend K{t, a) by for cr > r. It is convenient to define 

E 



T > cr > 0. 



(3.40) 



J(t,(t) 



W 



+ W{a -to){t -To)- 



with this notation, (3.39) becomes 



(3.41) 



(t - cr)J(T, cr) 

Ri{T)Ri{a) ' 



and we may also remark that 
(3.42) 



i?i(r) = ^J{t,t). 



Remark 3.1. The function K{T,a) > if r ^ 1: indeed the only possibility for 
K to be strictly negative is when the product (cr — to)(t — tq) is strictly negative, 
and cr < r; therefore, a is smaller than tq and r is larger than tq. Relations (3.25) 
and ( 3.33 ) imply that W^{to — a){T — tq) is estimated by Cr. Therefore, if r <C 1, 
if (r, cr) is nonnegative. 

3.4. Validation of the earliest asymptotic. With the help of the kernel K, we 
consider now the following problem: to find an interval [0, ri] and a mapping Si 
from this interval to R such that Ri + Si solves ( 3.17 ), with the initial conditions 
( p^ ) and ( |3l9| ). As Ri{0) = R{0) and 7?i(0) = i?(0), 5i(0) and 5i(0) have to 
vanish; therefore (3.17) can be rewritten as the following integral equation 

(3.43) Si{t)=Ci{Si)(t)+Gi{t), 

where Gi is the function 

K{T,a){2aRi + Ri){a) da, 

E{1 - ef 



(3.44) Gi(t) = 
and Ci is an integral operator defined by 

Ci{Si){r) ^ r {K{r,af^' 

(3.45) ^0 



R\ 



E 



(Ri + Si) 



Si 



(^) 



2a 



dK{T, cr) 



da 



da. 



Define 
(3.46) 



Ti = 77' 



1 < 71 < 2. 



Our purpose now is to prove that (3.43) has a unique solution on [0, ri] thanks 
to the strict contraction principle. 

We equip the space of continuous functions on [0, ti] with the norm 

(3.47) IIS'ill = max{|S'i(o-)|/i?i(o-) : < o" < n}. 

The choice of the weight Ri in the norm is natural since we expect Ri to be the 
principal part of the solution; thus we expect that the relative error (R — Ri)/Ri 
will be small: our norm measures precisely this relative error. 
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In order to apply the strict contraction principle, we estimate certain functions 
through a sequence of technical calculations. 



Lemma 3.2. For all large enough k, the expression 

1 r K{T,a) 

is bounded on [0, 1]. The bound will be called henceforth S. 



(3.48) 



/(r) 



■ da. 



Remark 3.3. The expression /(r) controls the nonlinear term in the integral equa- 
tion ( 3.43 ), whose detail is given in ( 3.45 ). 



Proof. We use the explicit expression of K to perform an estimate of /: 



/(r) 



1 



Rlir) J, 
We introduce the notation 



E 
W 



(3.49) 

and the change of variable 
(3.50) r = To 

The integral /(t) is now given by 

1 



Rtia) 

K = Ve/w, 



+ Wic7 -To){t -To) 



da. 



To 



KX. 



I{r) 



J{t,t) 
1 1 



^ {T-a)J{T,a) 



da 



-/(o-^o-)^ 



El + y 

We use the obvious inequalities 

(3.51) \l + xy\<Vl 



Tq/k 



(1 



X 



2^2 



dx. 



(3.52) 

to infer that 



\nT)\ < 



1 



\y\ 



\y-x\ < \x\ + \y\, 



dx 



1 



1 



(1+^2)3/2 E^T^J_,^,/,[l+X^fl^ 



\x\ dx 



It is now clear that |/(t)| is bounded independently of y, i.e. of r by a certain 
number S. □ 

We shall prove now that for large enough k, we can apply the strict contraction 
principle to the mapping 5*1 i-^ £i(S'i)-|-Gi on th e interval [0, n], defined by ( 3.46 ); 
for this purpose, we use the norm defined at (3.47), and we show the following result 



Theorem 3.4. For all 71 e (1, 2), and for all small enough p g (0, 1) there exists 
ko > such that for all k>ko, the mapping Si 1-^ Ci{Si) + Gi leaves invariant the 
ball of center and radius p and is a .strict contraction in that ball. In particular, 
if R is the solution of (3.17), (3.18) and (3.19), we have the estimate 

(3.53) ^ - ^M-^^M 



max 

0<T<77~'1 



< 



Rlir) 

The proof of this result depends on several estimates given in successive lemmas. 
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Lemma 3.5. For all 71 G (1, 2), there exists a constant C such that 
(3.54) IIGill < ^(774^1-4 + 77^1). 

Proof. By an integration by parts, 

Gi = 2aK{T, 0)i?i(0) + ^ ^2a^(T, cr) - if(T, ct)^ Ri{a) da. 



We estimate |Gi(T)|/i?i(r): we first observe that 



K{t,0)Ri{0) 



E 
W 



WTfiij - To) 



Ri{r) Riir) 

We estimate R\{t) fr om b elow hyEjW; we also observe that W'^tq is bounded, 
thanks to estimates ( 3.25 ) and ( ^3.33| ); therefore we have the following estimate, 
where we have used again remark 3.1: 

K{t, 0)i?i(0) 



(3.55) 



< r(l + VFVot/^;) 0(7/'^!). 



i?i(r) 

Next step is to calculate dK{T,a)/da: we use formulas ( ^.39| ), ( 3.40 ) and ( 3.41 ) 
and we find that 

dK 

— (r, cr) = £i(t, cr) +L2{t,cj) +L3(t, cr), 
acr 

where the Lj's are respectively given by 



Li{T,a) = 
L2{T,a) = 

i3(T,cr) = 



J(r,a) 



i?i(r)i?i(a)' 
(T-a)(a-To)J(T, (t)W^ 
i?i(T)i?3(a) 
(r-a)(r-ro)l¥ 



Lj{T, a) \ i?i(cr) dcr. 



i?i(T)i?i(a) 
Our aim is now to estimate the expressions 

The first expression /i(t) is rewritten with the help of the change of variable (3.50) 
and becomes 

ry 

h{r) = 

which we estimate as follows; 

h{r) < 



1 + J-r„/n 



\l + xy\ dx, 



y\ + xP' dx. 



TO Ik 



Since tq/k < 1 for fc large enough, we can see that 

inin(y ,0) 



—= / V^ + x^dx <V2 

^/l + y'^ J-to/k J-to/k 

y2p + 

V K 



dx 



i(y,o) Vl + 



dx 



< 



Therefore 



/i(t) <V2t. 
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With the change of variable ( 3.50| ), we may write the expression l2iT) as 

{y — x)\x\\l + xy\ dx 

(1+2/') J-ro/. 



h{r) 



1 + x^ 



We use ( ^.51 ) again together with 



< 1, 



and we find 



hir) < 



V^+y^ J -to/, 

But for k large enough and r > 0, 



{y — x) dx < 



2VE^l+{T-To)yK^ 



Vl + (t - To)2/k2 > — ; 

indeed, if t > 2tq, t — tq > t/2, and the inequality is clear; on the other hand, 
if T < 2to, for k large enough tq < k, and the inequality also follows. Therefore, 
there exists a number C such that for all large enough k and all r in [0, Ty'^'^J the 
following inequality holds: 

The third expression is handled as follows: 

1 r {t - a)\T - To\W da ^ \t~to\t^W 
Ri{t)Jo Ri{t) 2J{t,t) ■ 

If < r < Tq + At, we use the inequality J(t, t) > E/W and we obtain 



h{r) < 



2E 



since |t — ro| < max(K, tq) = k for k large enough. On the other hand, for r > tq + k 
and for k large enough 

r 

<2, 



T - To 

and therefore, using the inequality J{t,t) > W\t — tqP, we obtain 

, ^ \T-ro\r^W ^ 

l3[T) < 7T7TT7 < T. 



2t-to^W 



Thus, we have shown that 
1 



(3.56) 



Riir) Jo 
There remains to estimate 



dK 



(3.57) 



We rewrite (3.57) as 
1 



Rlir) 



R 



Rt{a) 



Ri{a) da = 0(7^''^). 

^ f K{T,a)Ri{a)da. 
liV Jo 



Rt(a) 



E 

— +W{a^To){T~To) 



da 



14 



LjGTITIA PAOLI and MICHELLE SCHATZMAN 



and we find that thanks to remark p.l 
R 



(3.58) 



^ / \K{T,a)\Ri{a)da 
il'T) Jo 

-1- r K{T,a)R,{a) < ^ max(i?t(0), i?f(r)) = 0{rj*''^-'). 
^iv) Jo 

Summarizing (3.58) with ( 3.55| ) and ( [3.56[ ), we find estimate ( |3.54| ). 
Next lemma enables us to estimate ||£i (5*1)11 when ||5i|| < p < 1. 



□ 



Lemma 3.6. Assume \\Si\\ < p < 1- For all 71 e (1,2), there exists ko and C 
such that for all k > kg the following estimate holds: 



^i(^i) < 



2eEd 



(l_p)5 (l_p)3 

Proof. The easiest part is the estimate on 



^ (^-K{T,a)Si{a) + 2a^{T,a)Si{<j)^ da. 



We can see that the absolute value of this expression is estimated by 

dK , 



\Si\ 



2a 



da 



[r, a 



K{T,a) Ri{a)da. 



We recognize expressions which have already been estimated in (3.56) and (3. 
Therefore, it is immediate that 



(3.59) 



1 



dK 

-^K{T,a)Si{a) + 2a—{T,a)Si{a) ) da 



Riir) Jo 
< C|| 5i II (774^1-4 + r;^i). 

Next comes the slightly more complicated expression 

1 



(3.60) 



Ri{r) Jo 



K (t, a) tt; „ , „ , , da. 



(i?i + Sif{a) 

For all fc > 0, 2 e — is at most equal to 2£\ therefore, if ||5i|| <p< 1, then the 
absolute value of ( 3.60 ) is estimated by 

K{T,a) 



2eE 

RA^Jo Rl{a){l~pr 
Therefore, we can see that 
1 



da. 



(3.61) 



Riir) Jo 



, {2e~e^)E ^ 
K[T,a) — „ , , da 



< 



2eE5 
(T^ 



{Ri+Sif{a) 

The last and most complicated term contains the expression 

E E SESi \ 

XR,TS^^M^~Rf) ^"'^ 
which can be rewritten thanks to Taylor's formula with integral remainder as 

Sf{a) 



12E 



{Ri{a) + sSiia)) 



-(1 - s) ds. 
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The absolute value of this expression is estimated by 



Therefore, we find that 
(3.62) 



Riir)l ^^^^'"^ {{Ri+s.r m 



(l-p)5i?i(a)3- 

E 3ESi 



(a) da 



< 



GESp"^ 



□ 



End of the proof of Theorem \3.4 - The ball of radius p about will be invariant by 
the mapping S'l i-^ Ci{Si) + Gi provided that 

QE6p^ 

(3.63) 



2eE5 



{l-pf i^-pf 

Choose p small enough for the following inequality to hold 



QE5p 1 
5 ^ 2' 



Choose then k so large that ( |3.63| ) holds. 

Let us prove now that for an adequate choice of p and feg and for all k > kg, the 
mapping £i is a strict contraction: the easiest part of the estimate pertains to 



1 



i?i(r) 



dK 

2aiSi ~ Si)ia) — iT,a) - (Si - Si)ia)KiT,a) ] da 



and it is clear from the proof of estimate (3.59) that 

dK 



(3.64) Ri{t) 

< C\\Si - Si\\{7j'^^'-'^ 



2aiSi ~ Si){a) — {T,a) ^ (Si ~ Si){a)KiT,a) da 
oa ' 



V 



The second easiest term involves the difference 



2e- 



2£- 



(i?l + ^i)3 (i?i + ^i)3 

which we estimate thanks to Taylor's formula: 

\iSi-Si){a)\ 



6e 



ds < 



6e||5i -5i 



/o (i?i(a) + Siia) + s{Siia) - Si{aW " " P)'-Ri(^)' ' 
Therefore, the corresponding term in £i(5i) — £i(5i) contributes an estimate 

6sE6\\S,-S4 

The last and most complicated term involves the expression 
3Si 1 3Si 1 
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We use a Taylor expansion twice to rewrite this expression as 

£ j-^ _{Si + s{Si~ Si))ds'd; 



12(^1 - S,] 



Ci{Si) is estimated 



{R,+s'{Si + s{Si-Si))y 

Therefore, the corresponding term is estimated by 

12p||^i-^i|| 
(l-p)5i?i(a)3- 

Thus, the norm of the corresponding contribution in Ci(Si) 

by 

(3.66) 12pESl\S. S^W^ 

Therefore, if we summarize the estimates ( 3.64 ), ( ^3.65 ) and ( t3.66| ), we find that 
on a ball of radius p < 1 about 0, the Lipschitz constant of Ci is estimated by 

QeE5 UpES 



(3.67) 



ri7l 1 



' ' (l-p)4 (l-p)5- 

If we choose p small enough for 12pE6/{l — p)^ to be less than or equal to 1/2, it 
is clear that we can choose fco large enough for the sum of the remaining terms in 
(3.67) to be less than or equal to 1/4. 

Together with the conditions found above for the invariance of the ball of radius 



p about 0, we have shown the first part of theorem 3.4. We also infer from this 
proof that 

11^1 II <C(r;4^^-4 



(3.68) 



V 



71 



I ts la st assertion is an immediate consequence of the equivalence of ( 3.43 ) with 
( 3.17), the definition of the norm, and the fact that the initial data coincide. □ 



Remark 3.7. Let us observe that in the end of the proof of Theorem 3.4, we can 
take p arbitrarily small provided that fco is chosen large enough. 

We conclude this section by the 

Proposition 3.8. Let ji belong to (1,2), and let ti — rf'^ . Then the following 
equivalences hold: 



(3.69) 
(3.70) 



R{t) - Riij) over [0,ri], 
R{t) ^ Riir) over [?7^,ti]. 



Proof. The first statement is an almost immediate consequence of (3.68): we have 



i?i(r) 



= 1 



i?i(r)' 



so that on [0, n] 



Rir) 



i?i(r) 



- 1 



and (3.69) follows 
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In order to compare R and we write the differential equations that they 
satisfy: 



R + 2aR = -R 
E 



i?3 



i?i + 2aRi = -^3 + 2ai?i. 
Ri 

Therefore, if we subtract the second of these equations from the first, we deduce 
that 



R - Ri + 2a{R ~ ill) 



Ei\~ef E 
R + - — - 2ai?i. 



i?3 Rl 
The initial data vanish. 

We integrate because we want to estimate R — R\: 

R{t)-Ri{t) 



(3.71) 



exp 



(-2a(T-a)){-i? + 



E{l-ef E 



— — 2ai?i >(ct) da. 



i?3 Rl 

The next step is to estimate the integral on the right hand side of ( |3.7l| ). We 
decompose this integral into three terms: 



Ii — ~2a I exp(— 2a(T — (j))Ri[(t) da, 
Jo 

h = -i exp{~2a{T - a))R{a)da, 
Jo 



l3 = E f cxp(-2a(r-a)) 
Jo 



-R^ 



The first two integrals are very easy to estimate: for Ii an integration by parts 
gives 

h = -2aRi{T) + 2ae-2"^i?i(0) - Aa^ [ exp(-2a(T - a))Ri{a) da. 

Jo 

Thanks to ( |337| ), on [0,ri], i?i(r) = 0{t]^^-^). Therefore 

(3.72) =0(7?^^-i). 

For I2, the situation is even simpler since it can be readily seen that 

(3.73) \h\<{l+p) I Riia)cxp{-2a{T-a))da^O{r]^^'-^), 

Jo 

where p — \\Si\\. 

There remains to estimate I3 ; a straightforward calculation shows that 



i?3 Rl 



< 



here p = ||5i|| is estimated at ( 3.68| ), so that 
(l-e)2 1 



(3.74) 
Therefore 



2£ + p(3 + 3p + p^) 1 

(l-p)3 Rr 



i?3 Rf 
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where J is defined as 



But J can be more conveniently rewritten as 



J = j^ expHa(r-.))-^da. 

TOniently rewritten as 

J= exp(— 2q:(t — (T))i?i(cr) dcr, 
"'0 



which we integrate by parts. We find that 

J = i?i(T) - e-^"^Ri{0) + 2a(i?i(T) - e-2"^i?i(0)) 

+ 4a^ [ exp{-2a{T - a)) Ri{a) da. 
Jo 

We can see now that 

J = 7?i(r)+ 0(77^1-1). 

Therefore 

Mr) , 0(77^1-1 ' 



i?l(T) \ i?l(T) 



Since i?i is nonnegative, for r G [77'^, ti], Ri{t) can be estimated from below by 
Ri (ry'^) which is equal to 

W{r^^ - To) 



we infer from relations ( 3.33 ) and ( 3.34 ) that 

lim Ri(t]^) > 0, 

which concludes the proof. □ 

4. The case < 7r/2. 

We prove here the first theorem which justifies Moreau's rule for 6 < tt/2. 

Theorem 4.1. If < 7r/2, the representative point of the system enters region TZ^ 
at a time i — to + f/Vk, where 



(4.1) fr^ta.ii9VE/W; 

moreover, in the coordinates defined by the axes yi and y2 (see Fig. we have 
the following asymptotics for all t > t: 

C cxpUiit - t)y/k) 

(4.2) 0<yi(t)< ^^'^ ^ 

(4.3) y2(t) ~ (i"t)s(O)cos0. 



Proof. The differential equation satisfied by O defined by ( 3.28 ) is deduced from 
(3.S) and is given by 

(4.4) ^ = ^^^^^ ®(0) = 0- 
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Recall that the principal part 0i is defined by ( 3.29| ) and (B.30). Let us estimate 
Ti = — Oi: Ti satisfies the differential equation 

~R2 



Therefore, if we let p = ||S'i|| , and if we denote 



(4.5) 

we find that 



Ti 



</3ei 



Hence, for all r G [0, n], 

(4.6) (1 - miir) < e(T) < (1 + /3)ei(r). 

According to ( |3.68 ) and the definition of e, there exists ko such that for all 
k > ko: 



(4.7) 



/?< 1, 



< 



nil - d) 



Let T_|_ and t_ be defined by the relations 



1-/3' ' 1 + /3 

Thanks to condition ( [4.7[ ) and formula ( 3.32| ) r_|_ and t_ are well defined, and are 
given by 

^ tan(^/(l - 13)) - WtqI^/E 



W 1 + {Wto/VE) tan(^/(l - /?)) ' 
E tan{e/{l + (3)) -Wto/Ve 



W l + {WTo/VE)t&n{e/{l + /3))' 

Therefore, as k tends to infinity, both t_ and t+ are equivalent to (^^/E/W) tan 9. 

The function Q is strictly increasing with respect to time; thanks to inequality 
(^^), there is a unique f e [t_, r+J such that 8(f) = 6. 

We know an equivalent of r_ and t+ as k tends to infinity: 

tan^ = 0{r]^). 



E 



(4.8) 



Together with (3.37), the above relation implies 

f(0)?7 



i?i(f) 



2 cos 



1 



and from (3.3S) that 
(4.9) 



i?i(f) 



s(0) sin( 



Proposition 3.8 implies R{f) ~ Ri{t) and i?(f) ~ i?i(f). 

We change coordinates now, taking the axis y2 along the second side of the 
convex cone K and the axis yi perpendicular to y2, and going out of K. The new 
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time variable is a translation of the natural time, denoted t' , and we set its origin at 
the time when the representative point enters region TZ^- We also let t = io + 
With these conventions, 

y(0) = r(f), 2/(0) = r(i) + ir{t)9{t). 

We use now the equivalents obtained previously: 

(4.10) yi(0) = 0(r,VVfc), 2/2(0) = 0, 

(4.11) yi(O)-s(O)sin0, 2/2(0) - s(0) cos ^. 

The second component 2/2 of y satisfies the ordinary differential equation 

2/2 = 0, 

so that 

(4.12) y2{t') ^t's{Q) cos e. 

The first component 2/1 of y satisfies the following ordinary differential equation 

(4.13) yi + 2aVkyi + kyi = 0, 



as long as j/i > 0. The explicit solution of (ITS) with initial data ( 4.1C ) and ( 4.11 ) 
is given by 

y^it') = 2/1(0)^ , + yM^ ; ^ . 

Since 2/i(0) is non negative, yi{t') stays non negative for all t' > and we have the 
following estimate on the first component of y: 

□ 



(4.14) 0<2/i(O< 



Thus, we obtain the conclusion of this section as the following Theorem: 

Theorem 4.2. Let Hi be the orthogonal projection on {xi — 0}, and let II2 be the 
orthogonal projection on {xi cos^+X2 sin^ = 0}; see Fig. As k tends to infinity, 
Uk converges uniformly on compact sets of to Uoo given by 



UoQ [t) 



u(0) + ffliu(0) tfO<t<to, 

{t - to)Il2lIiu{to) lfto<t. 



Proof. The initial part of the motion is described thanks to (2.6) and (2.7). Es- 
timate (4.14) proves that yi{t) tends to uniformly on compact sets of ]to,oo); 
relation ( 4.12 ) enables us to conclude. □ 

5. The second asymptotics 

In the case > 7r/2, we need a new asymptotic, and an estimate which is based 
essentially on the use of Lyapunov functionals, and which will be proved in Section 

i . . 

We restrict the choice of the exponent 71 in the definition of ti by assuming that 

(5T) 
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The reason for this choice is the following: if (5.1) holds, then the term E{1 — 
e)^ / R{Tif' is of order rf^'^'<^ which is small relatively to R{ti), according to the 
analysis of |]^: indee d, t he following equivalents of R{ti) and R{ti) are a conse- 
quence of proposition |3.8| : 

(5.2) i?(Ti) ~ s(0)r;Ti-\ i?(Ti)~ 5(0)77-1. 
Let ^ be such that 

(5.3) 0<C< 1/161- 
We define the time T3 by 

(5.4) T3 = C 111(1/77). 

We use the notation rs, because we will define below an intermediate time T2 

between ti and T3. 

We claim that the solution of ( 3.17 ) on the interval [ti,T3] is very close to the 
solution of 

i?2 + 2ai?2 + i?2 = 0, i?2(ri) = i?(Ti), i?2(Ti) = i?(ri). 
Let us define two kernels K2 and H2 on IR+ by 

K2{t) 



2VA ' 



2^/A 

We extend K2 and H2 to by 0. Therefore, R2 is given explicitly for r > ti by 

(5.5) i?2(T) = K2iT ~ Tl)i?(Tl) + i/2(r - Tl)il'(Tl). 



In order to substantiate our claim, we argue as for theorem 3.4: write R = 
^2 + '52 ; then 5*2 is a solution of the integral equation 

It is convenient to denote 

£2(^2)= r K2ir-a) J^^^^'l da, 
Jri (R2 + S2y{(y) 

whenever R2 + S2 does not vanish over [ti, T3]. 

Let us prove that R2 never vanishes over [ri,+oo): thanks to the inequalities 
> 6 > S,2, the functions K2 and H2 are positive for r > 0, -^2(0) is equal to 1; 
R{ti) and R{ti) are strictly positive. Thus the positivity of R2 is clear. 

On the space C'^([ti, ra]), we introduce the norm 

(5.7) 11^211 =sup{|^2(T)|/i?2(r) :Te [Ti,r3]}. 

We remark that £2 is well defined on the open ball of radius 1 about in the 
norm (5.7). 

We prove that C2 is a contraction on an appropriate ball, which will lead us to 
validated asymptotics for R on the interval [Ti,r3]. 

Theorem 5.1. For all p G (0, 1), there exists fci > such that for all k > ki, C2 

is a contraction from the ball of radius p (relatively to \\ \\) about to itself. 
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Proof. We will show in Lemma 5.2 that the expression 



i?2(T) J,^ Rlia 

tends to as /c tends to infinity, uniformly on [ti,T3]. 
If ||S'2|| <P < 1, then 

and in consequence, 

\\C2iS2)\\< ^,^^~'l' sup /(r). 

Let us estimate ||£2('S'2) — '^2(52)11 when max(||S'2||, ||S'2||) is at most equal to 
p < 1. We can see that 



2 



||/:2(52) - £2(52)11 <^§^—^ sup Iir)\\S-2-S2\\. 

Therefore, for k large enough, £2 is a strict contraction from the ball of radius p 
about to itself. □ 

Let us prove now the estimate announced on /: 
Lemma 5.2. The following estimate holds for I{t) on the interval [ti,T3\: 
(5.8) /(T) = 0(r;4-37i+^4(i+«i)). 



Proof. The integral / is analogous to the one defined in ( p.4§| ). 
We define T2 by 

(5.9) 

and we consider three cases: 

• Ti < T < 2ti: in this case H2 cannot be neglected relatively to K2. 

• 2ri < T < T2: in this interval, the dominant term in R2 will be R{ti) K2{t ~ti) 
and an elementary computation shows that this expression vanishes for r — 
Ti + (ln(^2/'^i))/(Ci ^ C2). Thus i?2 crosses approximately at a time r2/2. 

• T2 < r < T3: the last leg of the journey, since K2 is dominant and in K2, the 
term involving exp(^i(T — ri)) is dominant. 

Before proving these estimates, we observe that there exist positive numbers M 
and m such that 

(5.10) Vr e IR+, A'2(t) < Mr, 

(5.11) Vre[0,r2], K2{t) > mr. 

We tackle now the three separate sub-cases in detail. 
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5.1. First interval: re [ri,2Ti]. We remark that 

(5.12) R2{r)>R{n)H2{T-T,). 
Therefore, we can estimate /(t) as follows: 

(5.13) 0</(t)< r rfa. 

i?2(T)-R3(ri) Hi{a-Ti) 

We observe that over [ri, 2ti], 

(5.14) i/2(r-Ti)-l = o(ri), 

and we use ( ^.10 ). These observations imply the following inequalities: 

1 r K2{T-a) ^ 

da 



i?2(T)i?3(ri) 7,^ H|(a-ri) 



< C . / (r - a) da 



< C- 



i?^(Tl) 

and thanks to (f3.2|), the definition ( 3.46| ) of ti and condition (5.1), we obtain 
(5.15) /(r) = 0(7^4-2^1^^ 

5.2. Second interval: t e [2ri,r2]. We cut the integral / into two pieces: one 
piece from ti to 2ti on which we work essentially as in the previous sub-case, and 
a piece from 2ti to T2 on which we work differently. More precisely, on [ri, 2ri], we 
observe that R2{t) > R{ti)H2{t — ti), and on [2ti,T2], R2{t) > R{ti)K2{t — ti). 
Therefore, 

1 r-K2ir-a)^^ 



i?2(r) J,^ Rlia) 
r2Ti 



(^^3^371^^^ 



- R2ir) 

We estimate 1^2 (t — ti) from below by arguing that K2 increases from to a 
maximum, and then decreases exponentially fast to 0. Therefore, for all small 
enough 77, there exists r( tending to infinity such that K2{t[) = K2(ti). Moreover 
on the interval [ri, t[] K2{t) is greater than or equal to K2{ti). Thus, for all large 
enough fc, if 2(1') > K2{ti) on the interval [ti,t{ — ri], and therefore 

Vr e [2ri, r(], R2{t) > R{ti)K2{t - n) > R{t^)K2{ti). 

Thus, we obtain thanks to (t5.2|) 



i?2(r) ^ RUa) 

For the other piece, we estimate R2{t) from below by R{ti)K2{t — ti), and we 
obtain 

1 r if2(T-o-) , 1 r K2{T-a) , 
5.17) — — / \, , ^ < ^ / '-da. 

^2 T Ari " i?2(T)i?3(ri)i2ri i^l '^-n 
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We use (5.10) and ( ^.11 ) in the integral term of the right hand side of inequahty 
( ^.17 ), and we infer that 



1 

i?2(T)i?3(ri) 72ri Kl{a-n) 

c 



(5.18) 



< 



< 



i?2(r)i?3(Ti) Ari (^-Tl)3 

C{t - n) 



da 



da 



< 0(77^-2^1). 



The combination of (5.16) and (5.18) yields 

(5.19) /(r) = 0(77^-371)^ 

5.3. Third interval: r € [t2,T3]. We cut now / into three pieces, relative to the 
subintervals [ri,2ri], [2ri,r2] and [12, r]. 

On the last two pieces, we observe that t is far from ri, and we use the estimate 
from below 

(5.20) i?2(r) > 7?(ri)i^2(r-ri). 
Moreover there exists C such that for r > T2 and k large enough 

(5.21) if2(T-ri) > Cexp(ei(r-ri)). 

On the first subinterval, i.e. a g [ri,2ri], we use inequality ( ^.10 ); relations 
( p^ ) and ( ^.14D imply that 

thanks to (5.20) we can see that 



Mr) 



2ri 



K2iT-a) 

Rlia) 



da < 



1 + o(ti) 



RiTl)3R{Tl)K2{T-Tl) Jr, 



r2Ti 

/ M(t - a) da. 



Thanks to (5.21) and the asymptotics (5.2), we obtain 



da < Cr^^-'-'^ exp(-a(r - r,))r,r. 



i?2(r) Jr, Rl{a) 
Since t < and exp(— ^1X3) = r/''^^, we get finally 



da < C7y4-27i+C6 1^(1/,^). 



i?2(r) Jr, Ri{a) 



Relations (5^) and (5J) imply that 

4 - 271 + C6 > 0. 

We observe that 

(5.22) Va > 2ri, i?2(a) > R{Ti)K2{a - n), 

and we use estimates ( ^.1C| ) and ( ^.11 ); therefore 

f"^^ da ^ Mt3 



R^{Ti)J2r, m^{a~Tif - 2i?3(ri)TO3r2 
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Now, thanks to ( ^.20D and (| .2lh, we obtain 



Mr) 



Ki{T-a) 



da < ^ 



■exp(-67-3) < C77''-2'>'i+«i ln(l/7y). 



/2n " R{Ti)^rf 

Let us consider the third piece: we use now estimate (^-21 ) on the denominator 
of integrand; since K2{t) < Cexp(^ir), and thanks to (5.22), we have 

'•^ K^ir ~ a) 



1 

< 



Rl{a) 



■ da 



C 



exp(^i(T - cr)) 



da 



i?2(T)7?3(Ti) Jr^ exp(3a(^ - Ti)) 

and we conclude that the following estimate holds: 

(5.23) I{t) = 0(77''+«i"27i ln(l/77) + r;4(i+«i)). 

We have to keep the two terms in the above expression, since we have no way to 
ascertain the order of the exponents of 77. 

Whe n we compare the exponents in ( 5.15 ) and ( 5.19 ), we fi nd th at the exponent 
of ? 7 in ( 5TT9I ) is the smaller; when we look at the exponents in ( 5.23 ) to the exponent 
in ( 5.19 ) we fin d th at 4 + — 271 is strictly larger than 4 — 371, and this leads to 
the conclusion (5_^). □ 



We state now the main result of this section: 
Proposition 5.3. The following estimates hold: 

(5.24) R{r) ^ R2{t) uniformly over [tijTs], 

(5.25) i?(r3) 



- '/ ; 



(5.26) 



2^/A 



2VA 



Proof. Theorem 5.1 implies the uniform equivalence (5.24), and (5.25) is an imme- 
diate consequence of (5.24). 

Let us prove an estimate of the derivative i? at T3: 

^(0) 



i?(T3) ~ i?2(r3) 



(5.27) 

We observe that 

R{r3) = R2{r3 

Therefore, 

|i?(r3)-i?2(r3)| <C 



2VA 



'{i+ai 



dK2 
dr 



[T3 - a) 



{R2 + S2){ar 



da. 



dK2 



dr 



(t3 - cr) 



Rl{a) 



da. 



There exists a constant C such that for all fi > 

dK2, 



dr 



ia) 



< Cexp(^icr). 
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We use the method which gave estimate (5.23): we cut the integration interval into 
the three subintervals [ri,2ri], [2ri,r2] and [r2,r3], and on each of these subinter- 
vals, we estimate R2 from below exactly as in this calculation. Details are left to 
the reader, and we obtain 



exp(^ig) 



The equivalent of i?2(T'3) is obt ained immediately fro m th e explicit formula ( p. 5]) 
for i?2 and the equivalents (|).2|). Hence we infer that (5.27) holds. □ 



6. The final asymptotics 

In this section, we show that for large enough times R{t) is bounded from above. 
In view of (2.11), this estimate will enable us to show that the angular velocity is 
bounded from below, and hence, the polar angle Q will cross through 9. 

Theorem 6.1. There exists a strictly positive number Rm and a time T4 such that 

Vt > T4, i?(r) < Rm. 

Proof. Denote 



With these notations, equation ( 3.17 ) can be written 
(6.1) x = Mx + N{x). 



We observe that in the domain (0, 00) x R, (3.1) has exactly one critical point at 

1'), R.^{E{l-eff\ 

This critical point is attractive, as an examination of the linearization of ( |6.l[ 
around Xc shows. Moreover, there is a Lyapunov functional given by 

(6.2) F{x)=xl + ^^^''^\ xl 



Therefore, given x{t) with xi(t) > 0, we can see that for all r' > r, F{x{t')) is 
at most equal to F{x{t))^ and in particular, x(t') remains bounded. We see that 
when T tends to infinity, x{t) tends to the critical point Xc. 
The spectrum of M is {Cij?2}; therefore, the matrix 



/■oo 

Q— exp(sAf*) exp(sAf) ds 
"'0 



is well defined, symmetric, positive and definite. In particular, if Ai is the smallest 
eigenvalue of Q and A2 is the largest eigenvalue of Q, 

(6.3) Xix*x < x*Qx < \2X*x. 
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If we let x{t) = e**'^a;o: we observe that 
d d 

—x{t)*Qx{t) = — exp(iM*) exp(sM*) exp(sM) exp(tM)a;o ds 

dt dt Jq 

d [°° 

~ 'dt J '^'^Pi^M*) e'xjp{sM)xQ ds 
= -x{tyx{t) < ~X2^x{t)*Qx{t). 

As xo e and < > are arbitrary in the above calculation, we have proved indeed 
that for all x e 

(6.4) 2x*QMx < -X^^x*Qx. 

Since x = Mx + N[x) we have the inequality 

^^*Q^<-^^^2(x*QxY'^^2\N{x)\. 
dr M 

We seek a number R such that if R{t) > R, then 



Indeed, in order to satisfy (x5), it suffices to have 

2^,\N{x)\{x*Qxy/'<^, 

or equivalently, 

{x*QxY/^ 



\N{x)\ < 



But |iV(a;)| = E{1 - ef/R^ and |a;| > R, so that, with the help of (|6j), it suffices 
to satisfy 

E \\''^R 



I.e. 



We shall show now that if we choose R such that 



AEX 



3/2' 



i? > max ^ j , i?, 

then there exists r4 such that 

(6.6) R{ta) = R. 



Indeed, we know from (5.25) that Rir^) ^ 1, and that the limit of R{t) as r tends 
to infinity is Re, therefore, R{t) must cross R. We denote by T4 the smallest time 
in [r3,oo) such that (6.6) holds. 



On the interval [t3,T4], the differential inequality ( |6.5| ) implies 
{x*Qx){Ti) < {x*Qx){T3)exp{~{Ti -r3)/2A2), 
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< ln(a;*gx)(r3) -ln(a;*ga;)(r4). 



whence 

2A2 

But {x*Qx){Ti) > XiR^, and we obtain the inequahty 

T4 < T3 + 2A2 [ln(a:*ga:) {T3) - ln{XiR^)] . 
In particular, there exists C such that 

r4<r3 + Chi(l/77). 

We also need an estimate on R{Ti). We first show that it is less than or equal to 
0. By ( 5.26| ) we know that R{t3) < 0. Denote by (tsjTs) the connected component 
of {t > T3 : R{t) < 0} whose boundary contains t^. 

If Ts = 00, it is clear that R{t4) < 0. Assume that r5 < 00 and that R{t4) > 0; 
then Ts < T4 and R{tc,) vanishes. 

We infer from differential equation ( p. 17 ) that 

Rir,) = -R{r,) + 

but -R(t5) > Rc, because i?(T5) > R{t4) = R; therefore 
(6.7) i?(T5) < 0. 

On the other hand, as R{t) is negative on (tsjTs) and vanishes at rs, a straightfor- 
ward sign argument shows that 

Rin) > 0, 



which contradicts (6.7). 
Now, we prove that 



This will be a consequence of the following inequality for all r > ri and for all large 
enough k: 

(6.8) R{t)-^,R{t)>0. 
We observe that 

When we integrate this differential relation, we find that 

(i?- Cii?)(T) = exp(e2(T - n)){R - Cli?)(Ti) 



For k large enough, the equivalences ( ^.2| ) show that (7? — ^ii?)(ri) is strictly 
positive, and (6.8) follows immediately. 

We infer now from ( |6.8[) and the sign condition on R{t4) that 

F{x{t4)) < F ^ R^ + E{1 - efR-^ + ilR^ . 
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Since the Lyapunov functional decreases along trajectories of the system, we obtain 
for all T > T4 the inequalities 

(6.9) < i?(r)2 < F, R{Tf < F. 

F 

This is the final estimate we needed before the conclusion. □ 
We can now state the following corollary relative to the existence of the time f: 
Corollary 6.2. There exists a time f G (0,cxd) such that &{f) = 9. 



Proof. We know from ( 2.11 ) that 8 is an increasing function of r; if there is a time 
f < T4 for which 8(f) = 9, the conclusion is clear. Assume otherwise; then, with 
the notations of (3.£), we can see that 

and the conclusion is also clear. □ 

7. The case 9 > 7r/2 

In this section we estimate from below the first time f at which 8(f) = ^; we 
expect that f will be comparable to ri, but this is not correct. Recall the definition 
e; in this definition, the exponent of rj is 



4VA/|ei|; 

define a number r by 

(7.1) r = min(7i,4\/A/|ei|). 
Now, we can state the following theorem: 

Theorem 7.1. If 9 > tt/2, then for large enough k, f > ti; if 9 = Tr/2, then there 
exists a strictly positive number C such that 

(7.2) f > C?7'"^'^(^~''''^i). 



Proof. We argue as follows: assume r < ri; we recall estimate 

Vre [0,Ti], (l-/3)ei(r) <e(r) < (l + /?)ei(T), 

where (3 is given by ( [4.5[ ) and ||5i|| = p satisfies ( ^.681 ). The assumption (5.1) 
implies 471 — 4 < 71, and thus ( 3.6§| ) simplifies as 

p = 0{rf^ +e). 



The definition (7.1) of r implies that 

/3 = 0(r7'-). 

Moreover, relation (|3.32|) leads to 



8,(r)^--arctan ,+arctan^. 

This relation implies immediately that 

lim 8i(ti) = ^, 
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and therefore, thanks to (4.6) 

lim e(ri) - |. 

If ^ > 7r/2, the last relation implies immediately that for k large enough, f is at 
least equal to ri . 

Assume now that 9 — 7r/2; now, the situation is more delicate, since none of the 
inequalities established so far implies an estimate on f. If f > ri, we are done. 
Otherwise, we shall estimate f from below. Already, relation (4.6) implies 

- 2(1 

or in other words 



— — arctan ——, + arctan — ;= > 



2 W{f-To) - 2(1 + C77'-)' 

which implies 

\/E 

arctan—— < 0{rf), 

W(t - To) 

and therefore 

r-ro > Cry2-^ 

Thus, we have shown that 

f > Crf-''. 

If 7i > 2 — r, the relations 

Ti^rf^ >f> Cif-'' 
are contradictory for k large; therefore 

7i > 2 — r, k large ==> t > ri . 
Thus, we have shown ( |7.2| ). □ 



We deduce the following estimates from ( |7.2| ) and the asymptotics of sections ^ 
I and I 

(^^max(l-r,7i-l) < < ^ 

(7.3) „ 

\R{r)\ < - 
V 

We are able to show now the main result of this section: 

Theorem 7.2. For 9 > tt/2, as k tends to infinity, converges uniformly on the 
compact sets of to Uqo given by 



Uoc{t) 



M(o) + <niu(o) ifo<t<to, 

ifto<t, 
where Hi is the projection on the line {xi — 0}; see Fig. [J. 
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Proof. We go back to the original scales and time t — to + r/^/k; then 

Therefore, in coordinates yi,y2 (see Fig. |l|), we have the relations 

= 77i?(f)/Vfc, 
y2{t) = 0, 

= 77i?(f), 

We have also the estimate 

(7.5) i?(f) - eii?(T) > 0. 



(7.4) 



If f > Ti, (7.?;) is a consequence of (3.S). Otherwise, we observe that f belongs to 
[77"^, T i] for all large enough fc; therefore, we are able to use the equivalences ( 3.69 ) 
and ( 3.70 ), whence 

R{f) - ^iRif) ^ Ri{f) - ^Mf), 

which is valid because the dominant term in the right hand side of the above 
expression does not vanish; indeed, the expression ( 3.37 ) of i?i and ( 3.38 ) of Ri, 
we can see that 

Ri{f) - Cii?i(f ) - Vw~^iVWf - Cri-\ 



which implies (7.5); in the original coordinates, ( [7.5D translates as 

(7.6) yiit)-Vk^mit)>0- 



We infer from estimate (7^) that 

yiit)^0{l/Vk), 
yi(i)-0(l). 

In the coordinates yi and y2, the system ( |2.2D can be rewritten 
(7.7) yi + 2aVkyi + kyi^0 

as long as yi > and 

But the explicit solution of ( |7.7| ) with initial data ( |7.4| ) is given by 

exp(Ci(t - F)Vfc) - exp((6(t - 



y,{t)^yi{t)- 



■yi{t) 



2VAk 

a exp(e2(^ - t)Vk) - 6 exp(a(t - t)Vk) 
2VA 



If y(ti) is non negative, it is clear that j/i stays non negative for all time larger than 
i. If yi{i) is negative, we use ([7.6D: we estimate from below yi(t) by \/k^iyi(t), 
and after simplifications, we get 

yi (i ) > yi (t) %^ exp (ei {t -t)Vk). 



2VA 
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Therefore, (|7j) holds for all t > t 




yt > t, 



o{i/Vk), 



and 



Vi>t, y2{t) 



0{{t 



1 — max(l — r.71 — 1) 



which proves that in this case the limit of j/i and y2 is 0, as k tends to infinity. □ 
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